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I discuss recent developments in the study of cosmological limits on the 
Minimal Supersymmetric Standard Model (MSSM). In particular, I focus 
on the effect of neutralino-stau coannihilation on the relic abundance of 
neutralinos, and I give examples where the cosmologically derived limits on 
the supersymmetric parameters are relaxed, and one example (CP violating 
phases) where they are not. 

PACS numbers: 11.30.Pb, ll.30.Er, 95.35. +d 

1. Introduction 

The title of this talk is rather broad. Specifically what I will be talking 
about today is relic density constraints on supersymmetric models, and I 
will focus in particular on bounds on minimal Supergravity (mSUGRA). 
I'll begin with a brief reminder of the notation of super symmetry, then 
give an introduction to relic densities, how they are computed, and how we 
can use them to constrain models of low-energy super symmetry. I'll then 
discuss coannihilation in general and show how it dramatically relaxes the 
cosmological upper bound on the mSUGRA masses. Finally, I'll give an 
example of constraints (on CP violating phases in mSUGRA) which are not 
relaxed by the weakened mass limits. 



2. SUSY 

Since previous speakers have introduced the Supersymmetric Standard 
Model, I will give only a brief reminder, in order to present the particle 
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content and parameters and to set the notation. Recall that SUSY essen- 
tially doubles the particle content of the standard model (Table 1). Each 
fermion (in fact each fermion chiral state) has a spin-0 partner sfermion, 
the gauge bosons have spin-1/2 partner gauginos, and the degrees of free- 
dom of the Higgs sector, which now contains two Higgs SU(2) doublets, 
have spin-1/2 Higgsino partners. The four neutral gauginos and Higgsinos 
mix into "neutralino" states Xi, so that an arbitrary neutralino is a linear 
combination 



(1) 



where the B and are the partners of the U{1)y and neutral SU (2) gauge 
bosons and are linear combinations of the 7 and Z. The lightest of the the 
neutralinos xi tends to be the lightest supersymmetric particle, and in the 
models of interest for dark matter, xi tends to be S-like, i.e. ^ 1. 
Similarly, the charged gauginos and Higgsinos mix into two "charginos". 



SM SUSY 

fermions < > sfermions 

u,d,eL,eR. . . u,d,eL,eR. 



gauge bosons 
g,Z,j,W^... 

Higgs bosons 
Hi,H2 . . . 



gaugmos 

Higgsinos 
Hi,H2. . . 



Table 1. The SUSY partners of Standard Model particles. 



Along with the new particles of the MSSM come many new (soft SUSY 
breaking) parameters, including in principle separate mass parameters for 
all the sfermions, Higgs and gaugino mass parameters, and trilinear masses 
Ai of the Higgs-sfermion interaction terms, along with the supersymmetric 
Higgs mixing mass /x and its soft SUSY breaking counterpart B, and lastly 
the ratio of the two Higgs vacuum expectation values, tan/3 = V2/V1: 
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Sfermion masses 
Higgs masses: 
Gaugino masses: 
Trilinear scalar couplings: Ai hi Qi Uf H2 + . . . 
Higgs mixing masses: /j, Hi H2, B/j, Hi H2 

Higgs vev ratio V2/V1: tan (3 

If one includes flavor structure into the sfermion and trilinear masses, there 
are over 100 new parameters associated with softly broken supersymmetry[l]. 
This large number of parameters severely limits the predictive power of the 
the MSSM, and in practice, simplifications to the set of SUSY parameters 
are always made. One of the most popular and better-motivated choices is 
inspired by minimal Supergravity (mSUGRA) . In mSUGRA, several of the 
masses are taken degenerate, so that 

2 2 2 2 _ 2 /o^ 

'^h^ '^Jr " "''^i " "^^2 = ^0 (2) 

Ml = M2 = M3 = mi/2 (3) 

Ae = Ad = Au = ... = Ao (4) 

at the scale Mx where the gauge couplings unify. The parameters are then 
evolved to the electroweak scale using the Renormalization Group Equations 
to compute the low-energy spectrum. Due the RGE running, the masses 
of the sfermions depend on both niQ and my2. The parameters |^| and B 
are fixed by the conditions of correct electroweak symmetry breaking, which 
leaves 

"^0, mi/2,Ao,tan/3, sign(yu) 

as the free parameters of mSUGRA. I will subsequently concentrate on 
the mSUGRA model, but qualitatively similar results apply in the general 
MSSM. 



? ? 

JL JR 
o2 



m 



Hi 1 '^H2 

Mi,M2,M3 



3. The Neutralino Relic Density 

3.1. Relic Abundances 

The possibility of a significant relic abundance of neutralinos is partly 
due to R-parity, which is typically imposed on SUSY models in order to 
prevent rapid proton decay. Under this new symmetry, standard models 
particles have R-charge while their superpartners have R-charge — 1. 
Since R-parity is multiplicatively conserved, this implies that all vertices 
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must contain an even number of SUSY particles, and hence that the Lightest 
SUSY particle (LSP) is stable. Since it is stable over cosmological time 
scales, the LSP is a dark matter candidate, and we can use limits on its 
relic abundance to constrain SUSY models, as I'll describe next. 

Now, the very early universe was hot and dense; particles interacted 
rapidly, and the LSPs, which I will denote Xi were kept in chemical equilib- 
rium with the standard model particles in the thermal bath, primarily via 
processes in which two SUSY particles annihilate into standard model parti- 
cles, and the inverse processes in which standard model particles annihilate 
to produce two SUSY particles. 




Fig. 1. Keeping the x's in chemical equilibrium. 



At temperatures much larger than the mass of the the number density 
of x's was simply a spin factor times the number density of photons. As the 
universe expanded and cooled, the temperature eventually fell below the 
mass of the X; and the number density of x's began to drop exponentially. 



If this were the end of the story, it would be a rather dull tale: since the 
temperature of the universe today is about 3° K~ 2.5 x 10~^ eV, the number 
density of a, say, 100 GeV x would be suppressed vis-a-vis that of photons 
by a factor ~ exp{— 4 x 10^^}. I.e., there would be no x's left now. However, 
in an expanding universe, this conclusion does not hold, because at some 
point the x's fall out of chemical equilibrium with the bath. Specifically, 
this occurs when the x annihilation rate falls below the expansion rate of the 
universe, Pann ^ H. At this point, the x's cannot find each other in order 
to annihilate sufficiently fast for their number density to track the rapid 
exponential fall of (5). The number density of x's subsequently "freezes 
out" and simply falls with the volume of the universe, riy. ~ 1/V. 

To be more explicit, the number density of x's evolves according to the 
Boltzmann equation, 

dn 



= -3n^H - {a^nnv){nl - n^_^q ) (6) 
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Fig. 2. The x's fall out of chemical equilibrium. 



Here H = R/R, where R is the scale factor of the universe, and so the 
first term on the RHS of (6) simply represents the volume suppression of 
the number density of the x- The first term in parentheses describes the 
destruction of the x through annihilation, and the last term describes the 
production of x particles from interactions of the thermal bath. An approx- 
imate analytic solution to (6) is given by [3] 

^xh' - ^'7.. , (7) 



where = Px/Pc is the present mass density of x particles in units of the 
critical density pc required to close the universe, 5/ is the number of rela- 
tivistic degrees of freedom at freeze out, h is the current Hubble parameter 
H in units of 100 km/s/Mpc, and where the thermally averaged annihilation 
cross-section has been expanded in powers of (T/m^): 

(cTann^'> = a + b I — ] + . . . (8) 



The temperature at freeze-out is typically well below the mass of the Xi so 
that Xf = Tf/m^ ~ 1/20— 1/25, and eq. (8) is a good expansion^. Now, the 
key feature of (7) is that if the x annihilation cross-section is reduced, the 
X freeze out of chemical equilibrium earlier, when their density has had less 
time to track the exponential Boltzmann suppression (5), and the x relic 
abundance is larger. But a lower bound of 12 Gyr on the age of the universe 
(along with the assumption that Otot < 1) implies that i^toth'^ < 0.3. Or 



^ The temperature expansion is not good near s-channel resonances and just below 
important final state thresholds [2], where the cross-section can vary significantly with 
only small variations in the x energy. However, these occur only in a limited region 
of parameter space and are not significant for us here. 
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in other words, a lower bound on the age of the universe implies a lower 
bound on the x annihilation rate, and this is the feature we will primarily 
exploit to constrain SUSY models. 

3.2. Relic Density Constraints 

Let us now turn to the MSSM and see what these cosmological con- 
siderations do for us. The lightest supersymmetric particle is typically the 
lightest neutralino, and in many models, including mSUGRA, the lightest 
neutralino is a quite pure bino B. In the early universe, binos annihilate 
primarily via sfermion exchange into fermion pairs. 

X f 



X f 

Fig. 3. Sfermion exchange dominates x's annihilation. 



Now, if the mass of the sfermions is large, then B annihilation in the 
early universe is suppressed, and ^~h'^ is raised. Prom the last section, we 
see that the lower bound on the age of the universe impUcs an upper bound 
on the sfermion masses, and hence on both mSUGRA parameters thq and 
my 2. These limits are nicely complementary to those coming from direct 
searches for SUSY particles, which typically give lower bounds on the SUSY 
mass parameters. 

The cosmological limits can be translated into {my2jiT^o} plane [5], 
shown in Fig. 4. The light-shaded region corresponds to 0.1 < < 0.3; 

the area above this region is excluded. Below this region, so 
that another component of the dark matter would be required. This latter 
is not a bound in the same sense as the upper limit, since we don't know for 
certain that any of the dark matter is composed of neutralinos. In the nar- 
row chimney near my2 = llOGcV, rriy. ~ myj2^ and s-channel annihilation 
through the Higgs pole can bring the relic abundance of neutralinos below 
0.3, regardless of the sfermion masses. In the dark shaded region, the LSP 
is the right-handed stau, which is excluded by the very tight limits on the 
abundance of charged dark matter [3]. 

In Fig. 4 we also display current experimental limits: the light dotted 
contour represents the bound from searches for sleptons at LEP, while the 
the dashed line is a chargino isomass contour of 95 GeV, which approximates 
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Fig. 4. The light-shaded area is the cosmologically preferred region with 
0.1 < < 0.3. In the dark shaded regions in the bottom right of each panel, 

the LSP is the fa, leading to an unacceptable abundance of charged dark matter. 
Also shown are the isomass contours m^± — 95 GeV and rrih = 95 GeV, as well as 
an indication of the slepton bound from LEP. 



the LEP 189 chargino bounds at large mo- Note that the chargino bound 
excludes almost all of the Higgs pole chimney. The most significant experi- 
mental bound at this value of tan (3 comes from Higgs searches at LEP. The 
dot-dashed contour represents a light Higgs mass of 95 GeV, which approx- 
imates the Higgs limit from LEP189, and the bulk of the cosmologically 
allowed region is excluded. Now, the Higgs mass itself is sensitive to tan /3, 
and as tan /? is dropped, the dot-dashed contour moves quickly to the right. 
It is clear that for some value of tan /?, the Higgs contour moves to the right 
of the light-shaded region entirely, and this and lower values of tan /? are 
consequently excluded. The current bound at this value of tan P is around 
102 GeV [4], and these arguments imply a lower bound on tan /? of 3.7 (2.8) 
for fi < 0{fi > 0). We'll see in the next section that these constraints are 
weakened when we consider coannihilation. 



8 



4. Coannihilation 
4-1- The Basics 

So far, we have ignored interactions of the LSP with heavier SUSY par- 
ticles. Recall that the LSPs freeze out of chemical equilibrium when they're 
very cold {m^/T ~ 25), so that if the mass splitting between the LSP and 
the ncxt-to-lightcst super symmetric particle (NLSP) is 0(1), the number 
density of NLSPs at freeze-out is Boltzmann suppressed with respect to 
that of LSPs by a factor which is ~ exp{-25} < 10^1°. Therefore we don't 
have to worry about NLSP interactions. If, on the other hand, the LSP 
and NLSP are closely degenerate in mass, then the NLSP interactions near 
freeze-out may affect the LSP relic density. 

This produces two competing effects. First, the NLSPs freeze out of 
chemical equilibrium with the standard model bath at the same time as 
the LSPs and subsequently decay into LSPs, and so a significant NLSP 
abundance at freeze-out can increase the relic LSP density. Typically a 
larger effect is that since the NLSP interactions contribute to the exchange 
of particle number between SUSY and standard model particles (and can 
dominate, as we'll see below), the SUSY particles remain in chemical equi- 
librium with the thermal bath for longer and track the equilibrium down to 
lower temperatures, and this reduces the LSP relic abundance. 

How degenerate do the LSP and NLSP states have to be in order to 
produce a significant effect? Well, 

^NLSP . ^-Am/r . . ^-25(mNT,Rp/mT,sP-l\ ("g-) 

If the NLSP is 10% (5%) heavier than the LSP, this ratio is ~ ^(i). We see 
that unless the lightest states are highly degenerate, coannihilation will only 
be important if ctnlsp-lsp (or (Jnlsp-nlsp) > ctlsp-lsp- And mSUGRA 
(and much of the MSSM), they are!! Consider the temperature expansion of 
the thermally averaged cross-section (8). When the final state is a fermion 
pair (the dominant annihilation channel for a i?-like neutralino), a ~ mj. 
This dependence is due to the fact that one has identical Majorana fermions 
in the initial state [6] and is called "p-wave suppression". Since T/m^ is 
small at freeze-out, this suppresses the annihilation rate (and enhances the 
relic abundance) by an order of magnitude. Coannihilation cross-sections do 
not have such a suppression and are typically an order of magnitude larger, 
and the NLSP interactions can therefore dramatically reduce the SUSY relic 
abundance. These effects have been well studied in SUSY for Higgsino- 
like neutralinos [7], where there is typically a close degeneracy between the 
lightest and next-to-lightest neutralinos and the lightest chargino. What we 
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have found is that coannihilation is also an essential element in determining 
the cosmological upper bound on gaugino (B) like neutralinos, as well [8]. 

4-2. B — f Coannihilation 

Looking back at Fig. 4, we shouldn't be surprised that coannihilation 
may be important in mSUGRA. The upper bound on mi/2 occurs at the 
intersection of the ^^^^ = 0.3 contour with the the top of the region with 
mfj^ < rriy., i.e. at a point where the stau and neutralino are exactly de- 
generate! Generally in the MSSM ^, the cosmological upper bound on the 
mass of the B is saturated when the masses of the lightest sfermions are 
degenerate with m^. In mSUGRA, the three right-handed sleptons Trj/Ir 
and are the lightest sfermions and can all be close in mass to the LSP. 
We must therefore consider [2, 5] the effective annihilation cross-section 

where i,j = Tr, r*, Cr, e*,MR) Mr ^^"^ and where n = Yjf^i- The 
complete set of initial and final states contributing to (10) is given in Table 2. 

The dominant contributions to cJefr come from £r£r annihilation to gauge 
bosons, ^R^R annihilation to lepton pairs and £rX annihilation to a lepton 
plus a gauge boson. The final states with heavy Higgses turn out to be 
kinematically unavailable in the regions of interest. For further calculational 
details, see [8]. 



Table 2. Initial and Final States for Coannihilation: = r, e, /j.} 



Initial State 


Final States 


Ji Ji* 


77, ZZ,^Z, W+W-, Zh,^h,hh, ff, 




ZH, 7H, ZA, W+H", hH, hA, HH, HA, AA, H+H" 


Ji Ji 

'-R ^R 


t'ii 


•^R ^R ) ^ 7^ i 




4x 





^ The presence of s-channel heavy Higgs poles can provide a loophole when there is a 
small admixture of Higgsino in the lightest neutralino state. 
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In Fig. 5, we show the contributions to a = a + ^^Xf (see (7)) for 
TrTpJ annihilation. The top sohd contour is the total a for TrT^, while for 
comparison we display as a thick dotted line the equivalent total neutralino 
annihilation cross-section. As advertised, the stau cross-section is over an 
order of magnitude greater than that for the neutralinos, which is p-wave 
suppressed. Figures for — x and — annihilation show a similar 
enhancement over the x ~ X cross-section, and figures for other tan (3 and 
niQ are similar. 




Fig. 5. The separate contributions to the Tj^t* cross section a = a + i^bx for x — 
T/m^ = 1/23 and nio — 120 GeV, as a function of my 2- For comparison, the thick 
dotted line is the xx cross section. 

In Fig. 6, we display the contributions to o"eff as a function of the frac- 
tional mass difference AM = {ruf^ — m^/m^ between the neutralino and 
the stau. The thick solid contour shows the total agg, while for compari- 
son the thin solid contour gives the iTeff one would compute if one ignored 
coannihilations, i.e. a~~ + Here we've fixed mi/2 = 500 GeV and 

scanned upwards in mo, which increases AM. When the neutralino and 
stau are degenerate, the dominant contribution to comes from slepton 
annihilation. The ratio of the solid contours at this point is greater than 
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an order of magnitude, as above. As AM increases, the density of slep- 
tons becomes Boltzmann suppressed, and the slepton-slepton contribution 
falls with two powers of rf^ and drops below the slepton-neutralino contri- 
bution at AM ~ 0.07. This contribution in turn falls with one power of 
n^, and neutralino annihilation becomes dominant again at AM ~ 0.17. 
At large AM, the two solid contours and dot-dashed contour merge, and 
coannihilation can be neglected. Again, figures for other tan (3 are similar. 




0.1 0.2 

AM 



Fig. 6. The separate contributions to the cross section CTeff for x = Tjm^ — 1/23, 
as a function of AM = (rufj^ — m^)/m^, with (mi/2, tan/3) = (500 GeF, 3) 

Let's now look back at the {?7ii/2, ""t-o} plane and examine the effect of 
coannihilation on the cosmologically allowed region. In Fig. 7, we show 
the same area of parameter space as in Fig. 4, but now with coannihilation 
included. We see that the light-shaded area now bends away from the 
forbidden stau LSP region and creates a large allowed trunk which lies on 
top of the line ruf^ = m^. Eventually, for large enough my 2, the top of the 
trunk falls below the m^^ = line, but this doesn't happen until much 
larger values of my 2 and "^Oi as seen in Fig. 8. Some features of the new 
cosmologically allowed region to note: The upper bounds on mo and my2 
are relaxed to my2 < 1400 GeV and mo ^ 350 GeV, respectively. The 
width of the new allowed trunk is significant, from 30-50 GeV in mo for 
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mx/2 up to ~ 800 GeV. We've only shown the cosmologically interesting 
region for one value of tan /3, but the shape is very similar for all small to 
moderate tan /?. The position of the line m^^ = can vary somewhat with 
however, the width of the trunk above the line is quite insensitive to Aq, 
as is the upper bound on my2- This relaxes dramatically the cosmological 
upper bound on the neutralino mass in mSUGRA from about 200 GeV 
to close to 600 GeV. It will therefore take the reach of the LHC to probe 
the full cosmologically interesting region. Lastly, the bounds on tan f3 from 
combining the Higgs search limits with relic density constraints are now 
weakened, from 3.7 (2.8) to 2.8 (2.3) for /i < (/i > 0). 

200 



100 



a)tanp=3, |l< 




100 200 300 400 500 600 
^1/2 

Fig. 7. As in Fig. 4, but now including neutralino-slepton coannihilation. 



5. CP Violation and Electric Dipole Moments 

As an application, we'll now examine to what extent the relaxation of 
our upper limits on my2 affects constraints on CP violation in the MSSM. 
We'll start with a brief reminder of where CP violation arises in our model. 
Recall that in the MSSM, new CP violating phases 0^ and accompany 
the (in principle) complex parameters /x and Ai, introduced in the first part 
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2000 



Fig. 8. The same as Fig. 7, but extended to larger my2- 



of this talk. These phases then appear in the low energy Lagrangian in 
the neutralino and chargino mass matrices (in the case of 0^) and in the 
left-right sfermion mixing terms (both 6^ and 9a). The new sources for 
CP violation then contribute to the Electric Dipole Moments (EDMs) of 
standard model fermions, and the tight experimental constraints on the 
EDMs of the electron, neutron and mercury atom place severe limits on the 
sizes of 9^ and 9a [12, 13]. 

The EDMs generated by 6^ and 9a are sufficiently small if either 1) the 
phases are very small (< 10~^), or 2) the SUSY masses are very large 
{O (a few TeV)), or 3) There are large cancellations between different con- 
tributions to the EDMs. In mSUGRA, option 2) is forbidden by the relic 
density constraints, as we'll show next. Condition 3), large cancellations, 
does naturally occur in mSUGRA models over significant regions of param- 
eter space, including in the body of the cosmologically allowed region with 
"^1/2 = 0(100 — 400 GeV). These cancellations relax the constraints on the 
phases, but the limit on 0^ remains small, 6^ < vr/lO. 

To see why option 2) is cosmologically forbidden, recall that the SUSY 
phases contribute to the electron EDM, for example, via processes of the 
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following type: 




where selections and sneutrinos appear in the loop. These contributions 
diminish as the sfermion masses are increased, but this also shuts off neu- 
tralino annihilation in the early universe, which is dominated by sfermion 
exchange as in Fig. 3. The upper bound on ^^h^ then limits the extent to 
which one can turn off the electron EDMs by raising the sfermion masses. 
The combination of cosmological with EDM constraints in the MSSM and 
mSUGRA is discussed in detail in [14, 13]. 

To see the combined limits on 9^ and 6a in mSUGRA, we plot in the 
{^ixj^a} plane the minimum value of my2 required to bring the EDMs 
of both the electron and the mercury atom ^^^Hg below their respective 
experimental constraints (Fig. 9). These experiments currently provide the 
tightest bounds on the SUSY phases^. Here we've fixed tan/3 = 2, Aq = 300 
GeV and mo = 100 and scanned upwards in my2 until the experimental con- 
straints are satisfied. Due to cancellations, the EDMs are not monotonic 
in my2j however, there is still a minimum value of my2 which is allowed. 
Looking back at Fig. 4, we see that in the absence of coannihilations, there 
is an upper bound on my2 of about 450 GeV (though slightly smaller for 
this mo); an analogous figure to Fig. 7 for tan/3 = 2 shows that coannihi- 
lations increase the bound to about 600 GeV. Comparing with Fig. 9, we 
see that zone V is cosmologically forbidden, and that the effect of including 
coannihilations is to allow zone IV, which was formerly excluded. 

Note in particular that the overall upper bound on 6^ in this figure, 6*^ ^ 
0.047r, is not affected by coannihilations. This is because the largest 9^ occur 
in regions of cancellations, and these regions happen to lie at lower values 
of my 2, starting in zone II with my 2 < 300 GeV. Increasing m^j^^ from 450 
to 600 GeV is insufficient to bring the individual contributions to the EDMs 
to acceptable levels for the larger values of 9^ and significant cancellations 
are still necessary. Even taking m,y2 and mg at their maximal values from 
Fig. 7 is not sufficient to reduce the EDMs below their experimental limits, 
and so coannihilation does not affect the upper bound on 9^. 

^ The extraction of the neutron EDM from the SUSY parameter space is plagued by 
significant hadronic uncertainties [13], so that the inclusion of the neutron EDM con- 
straint does not improve the limits when the uncertainties in the calculated neutron 
EDM are taken into account 
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-0.05 -0.03 -0.01 0.01 0.03 0.05 



Fig. 9. Contours of rriy'^, the minimum mi/2 required to bring both the elec- 
tron and Hg EDMs below their respective experimental bounds, for tan/3 — 
2, mo = 130 GeV, and Aq = 300 GeV. The central light zone labeled "I" has 
m5^'2" < 200 GcV, while the zones labeled "11", "III", and "IV" correspond to 
200 GeV < mf/^ < 300 GeV, 300 GeV < m^^i^ < 450 GeV, 450 GeV < m'^f^^ < 
600 GeV and my™ > 600 GeV, respectively. Zone V is therefore cosmologically 
excluded. 



The bowing to the right of the contours in Fig. 9 is a result of cancel- 
lations between different contributions to the EDMs [14], and we can see 
that the effect is to relax the upper bound on 0^ by a factor of a few. As 
we increase ^40, the extent of the bowing increases, and larger values of 6^ 
can be accessed. This loophole to larger 0^ is limited by the diminishing 
size of the regions in which there are sufficient cancellations to satisfy the 
EDM constraints. In general, the regions of cancellation for the electron 
EDM are different than those for the ^^^Hg EDM, and the two regions do 
not always overlap. As 0^ is increased, the sizes of the regions of sufficient 
cancellations decrease; in Fig. 9, the width in of the combined allowed 
region near the 0^ upper bound is 40-80 GeV, which on a scale of 200-300 
GeV is reasonably broad. Larger permits larger 6^, but the region of 
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cancellations shrinks so that a careful adjustment of my2 becomes required 
to access the largest 9^. At the end of the day, values of much greater 
than about vr/lO cannot satisfy the EDM constraints without significant 
fine-tuning of the mass parameters. At larger values of tan /3, the upper 
bound decreases roughly as 1/ tan/3. See [13] for more details on the status 
of EDM and cosmological constraints on CP violating phases in mSUGRA. 

6. Summary 

In summary, constraints on the relic abundance of LSP ncutralinos place 
significant restrictions on the parameter space of the MSSM, and mSUGRA 
in particular. To accurately compute the cosmological upper limits on 
MSSM masses requires the inclusion of coannihilation effects, both for the 
case of a Higgsino and gaugino like neutralino. In particular, we have found 
that slepton-neutralino coannihilation greatly affects the neutralino relic 
abundance when the neutralino and slepton are closely degenerate in mass, 
as is the case in mSUGRA near where the cosmological upper bound on the 
neutralino mass is saturated. Including coannihilation effects significantly 
relaxes the cosmological bounds on my2^ "^o and m^, and reduces the com- 
bined Higgs + cosmology lower bound on tan (3, although the upper bounds 
on CP violating phases in mSUGRA are not relaxed. The reach of the LHC 
will be needed to be sensitive the full cosmologically allowed region. Lastly, 
although I did not discuss it here, similar effects are present for a gaugino 
like neutralino in the general MSSM. 
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